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1. Research. O b ject ive s 

Our research objectives, as stated in the work statement of our proposal, were: 

1. develop methods for providing low noise optical elements by interferometric 
means, including gratings, lenses, beam combiners, corrector plates, and other 
special elements, such as circular, eliptical and hyperbolic axicon lenses, including 
the following tasks: 

(a) apply the Lau effect in combination with spatial filtering to produce 
nonsinusoidal fringe distributions 

(b) compare the broad spectrum with the broad source techniques for HOE 
generation 

(c) seek methods for making HOEs in light that is simultaneously spatially and 
spectrally incoherent 

(d) using higher order analysis, study the aberrations of incoherent light HOEs, 
and methods to reduce aberration 

(e) determine methods for producing fringes under conditions of greatest 
possible incoherence, including determining any required departure from 
sphericity of the interfering waves 

(f) develop methods for producing HOEs having arbitrary fringe spacings, 
instead of only linear and quadratic, thereby producing aspheric HOEs, as 
well as arbitrary phase plates 

(g) determine the usefulness of sources with structured spatial or spectral 
distributions in HOE construction 

(h) assess the limitation on space bandwidth product of incoherent light HOEs 

(i) develop a comprehensive theory for the construction of HOE lenses in light of 
reduced temporal coherence 

2. Continue our basic investigations, with the intent of discovering new principles 
and techniques in the area of optical processing with light of reduce coherence. 

2 , Ac c omplishm ent s 

We have, during the past year, given attention to all of the above, with some of 
the objectives being investigated more intensively than others. Our principal 
accomplishments have been: 

1. We continued our work on the basic treatment of HOE generation with spatially 
incoherent light and developed a theory of broad source interferometry based on 
source distortion. We showed that for an interferometer to give broad source 
fringes, the interferometer must form distorted images of the source. The grating 
interferometer introduces exactly the correct distortion, the Mach Zehnder gives 
no distortion at all, hence has very limited broad source capability. 








2. We extended the theory of achromatic interferofmetry for HOE formation. We 
generalized the theory of the two zone plate interferometer for HOE formation. 

3. We developed a merit factor for describing how well an incoherent source can 
produce a HOE 

4. We developed a new method of spatial filtering using extended source incoherent 
light. The new system acts on the amplitude transmittance of the object 
transparency, and is therefore essentially a coherent system, even though the 
light is incoherent. 

5. We developed a new technique for exceeding the classical resolution limit in an 
imaging system (a process sometimes called superresolution). We showed that if 
we have a two-branch interferometer producing broad source fringes, and if we 
have an imaging system in each branch, and then, if we place an object 
transparency in one branch and reduce the aperture of this imaging system to a 
very small value while retaining a large aperture for the empty branch, then the 
resolution achieved is that corresponding to the larger aperture, even though 
light from the object passes through only the small aperture, not the large one. 










Introduction 

The construction of Holographic Optical Elements (HOEs) has recently 
become an important area ot research. Applications of HOEs such as heads up 
displays, supermarket scanners, and diode laser collimators have become 
increasingly widespread. The study of aberrations and noise of HOEs are two 
important topics of research. In this paper we discuss techniques for analysis and 
construction of low noise, zone plate HOEs. Low noise HOEs have been created by 
reducing the spatial coherence of the illuminating laser source 1 - 2 or by using a 
source of reduced temporal coherence 2 . This paper discusses the theory for 
analysis and construction of HOEs made by the latter technique, that is, HOEs 
made with a polychromatic point source. 

Optical systems have been designed to make holograms in polychromatic 
point source illumination 4 - 5 . The optical systems discussed by Katyl 4 result in either 
on-axis or off-axis holograms. The on-axis technique is limiting. We would like to 
generate holograms of arbitrary center spatial frequency. The off-axis techniqueof 
Katyl images a grating to the recording plane, resulting in a hologram that is noise 
limited by this grating. One would like to be able to produce off axis HOEs of arbitrary 
center spatial frequency, but not at the expense of imaging noise planes, such as a 
grating, to the recording plane. Or, stated somewhat differently, the plane of fringe 
formation should not be an image plane tor some physical structure, such as a 
grating, located upstream. 

The optical setup described by Collins 5 , employing a grating interferometer, 
was shown to be ideal for generating HOEs of this nature 2 . However, these zone 
plate HOEs were limited to interference between a plane wave and a spherical wave. 
Ideally, one would like to interfere waves of arbitrary curvature. This is a result of the 
aberrations inherent in a hologram when it is not used in exactly the same manner in 
which it was constructed 5 . Specifically, we would like to interfere a diverging wave 
with a converging wave. This is the only case that cannot be done cleanly with 
coherent light, since it necessitates at least one lens being introduced into one of the 
beams downstream from a pinhole, causing noise in the form of diffraction patterns 
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of dust or other scatterers on the lens. With the use of polychromatic illumination 
these diffraction patterns can be smeared into a negligible background. 

It is possible to construct zone plate lenses, as described above, with an 
arbitrary set of conjugate focal planes and arbitrary center spatial frequency, by 
generalizing the Katyl-Collins achromatic holographic system^. We shall see that by 
either moving or removing the final lens of Collins’ system (Fig. 1) we will be able to 
construct zone plate lenses, although of limited space bandwidth product, in 
polychromatic light. 



Fig 1. The system described by Collins, consisting of the grating interferometer and a four/two lens 
system illuminated by a plane wave. 

Etequ i re ments for Ach _rpma.tic ity 

The analysis for the construction of zone plate lenses in polychromatic point 
source illumination is formulated in terms of the zone theory of Fresnel in conjunction 
with the Newtonian lens equation. The radius of the m th zone of a Fresnel zone 
plate is given 

R m 2 = m: , F m = 1,2,3... (1). 

We see that the radius of the m th zone is directly proportional to the wavelength. As 
a consequence, the zone plate so constructed will exhibit a focal length F inversely 
proportional to wavelength, as in Fig. 2, resulting in wavelength dispersed point 
images. Reciprocity occurs when we construct a zone plate interferonletrically, by 








Fig. 2 A polychromatic plane wave hits an interferometrically recorded zone plate and is dispersec 
accorddmg to the formula F=F 0 (\ 0 /\). Where F 0 is the focal length at the midband wavelength \ 0 . 

interfering a plane and a spherical wave. If the spherical wave originates from a 
single point a distance F from the recording plane, the zones of constructive 
interference will have radii proportional to wavelength. However, if the point source 
is dispersed in such a way that the point source to recording plane distance is 
inversely proportional to wavelength, then each wavelength will make exactly the 
same zone plate. We seek a method by which we can make this wavelength 
dependence approximately constant over some wavelength band. From Eq. 1 we 
see that the only possibility for providing wavelength compensation is through F, the 
effective focal length of the zone plate being constructed. By manipulating the 
parameters that determine F, namely the radii of curvature of the two interfering 
beams, we can cause it to vary inversely with wavelength in an approximate fashion 
We know F is given in terms of the relation, 

1,'F = 1/r 1 + 1/r 2 (2), 

where r-| and are the radii of curvature of the two interfering spherical waves and 
may be functions of wavelength. We can generate arbitrary radii of curvature with the 
desired wavelength variahon by use of zone plate lenses in an optical system 

We obtain the radius of curvature of a wave by finding the location of the 



plane at which the source is imaged. This location is in general a function of 
wavelength, determined by the focal lengths of the zone plates involved in the 
optical system. Once this plane is located we back up or go forward from it to the 
appropriate location for the particular HOE design. The distance of this move is one 
of the radii of curvature described above, and is also, in general, a function of 
wavelength. This technique of analysis is best illustrated by example and will be in 
the following sections. 

One might notice that throughout this discussion no mention of actual 
physical pathlengths has ever been made. This comes about for two reasons. First, 
the grating interferometer, which splits and recombines the light, is inherently, 
minus all lenses, an achromatic system 7 . Second, both beams pass through the 
same optical elements in the system and because of this travel the same optical 
paths. The difference in sphericity of the two beams is generated by the use of zone 
plate lenses which allow a particular beam to be affected by the spherical phase 
transformation that is present on the plate, while the other beam is unaffected. This 
type of system, when analyzed using the above constraints, is used to produce zone 
plate lenses in polychromatic illumination. 

A Four Lens System 

Our analysis for the construction of HOEs with polychromatic illumination 
starts with the achromatic Fourier transforming system of Collins (Fig. 1 ). It was 
shown in Ref. 3 that off axis zone plate holograms could be made with this 
interferometer. These zone plates were limited to interference between a plane 
wave and a spherical wave. In this section we will analyze this interferometric system 
in a simple, very general fashion using the Newtonian lens equation for cascaded 
lens systems 8 . Then we will further generalize this system to produce zone plates 
with arbitrary conjugate focal planes. 

The interferometer of Fig. 1 consists of two beams, the object and reference 
beams, labeled O and R respectively. The object beam passes through four lenses. 
Lf L 4 While the reference beam passes through all four of these lenses, it is 
affected only by the two achromatic glass lenses L 2 and I 4 . This is a result of 
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selecting only the zero order from l_i and L 3 , which are zone plate lenses. Thus we 
have two lens systems to be analyzed under the requirements for achromaticity The 
easiest to analyze is the reference beam. The lenses l _2 and L 4 are in the afocal 
position so that when a plane wave enters from the left it exits as a plane wave to the 
right. Thus, the source, which is at infinity, is imaged to infinity. The radius of 
curvature of this beam we call r 2 and is in this case infinite. The object beam is more 
complicated. A point source at -infinity is imaged, by the four lens system, to a plane 
located by D as seen in Fig. 3. D is given in terms of the focal lengths of the four 
lenses by 

D = f 1 + f 2 -f 2 2/f 3 (3). 



D is a function of wavelength. This is a result of the two zone plate lenses Li and I 3 
whose focal lengths, fi and f 3 , vary with wavelength in the following manner 


f 1 = * 10 * 0 ^ 


(4) 








(5) 
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f 3 = *30*0^ 

Here the subscript zero is used to represent values of the focal length at the center 
wavelength Aq. Now we substitute these expressions into Eq. 3, allowing the 
wavelength, A, to be expressed as 

'\ = '\ 0 + AA (6). 

The source image location is now given by 

D = f 10 [1/(1 + AA/A 0 )] + f 2 - (f 2 2/f 30>< 1 + *a/a 0 ) ( 7 )- 

Using a narrowband approximation we write 

fl = fiot^'O + AA/Aq)] = fio (1 - AA/Aq) (8). 

Substituting Eq. 8 into Eq. 7 we obtain the final expression for D, 

D - (-2f 2 2 +f 2 f 30 )/ / f 3 o + (fi 0 + f 2 2/f 30)( 1 * ^A/A 0 ) (9) 

We note that the second half of the Eq. 9 is of the same form as Eq. 8, i.e. the 
second half is inversely proportional to wavelength. This is exactly the form for the 
radius of curvature which is required to achromatize the fringe pattern. Thus, r-|. the 
radius of curvature of the object beam is 

r i =(fio + f 2 2/f 30 )( 1 -AA/a 0 ) ( 10 ), 

where 

(1 - aa/a 0 ) - 1/(1 + aa/a 0 ) = Aq/A (11). 

With this approximate form the radius of curvature of the object beam satisfies the 
requirements for achromaticity. If a recording plate is placed in the plane D\ as seen 
in Fig. 1, where D' is just D minus r-| given by 

D' = (-2f 2 2 + f 2 f 30 )/f 30 (12), 

one is able to record an off axis zone plate in polychromatic illumination. Physically 
we can interpret this with the help of Eq. 13 which describes the resultant fringe 
pattern. The wave number, k, must be constant for all wavelengths at the plane 
located by D'. 






in 


2 + 2cos(kx 2 ) ( 13 ) 

The wavelength bandwidth of the illuminating source is limited since the analysis is 
only approximate. 

One can see from the above equations that Collins' achromatic Fourier 
transform analysis is a special case of the above analysis. If the radius of curvature, 
r-| (Eq. 11), is set equal to zero one obtains the condition 

f 10 = ’ f 2 2/f 30 ( 14 )- 

This is precisely the condition obtained by Collins. A radius of curvature of zero 
describes the source image plane as being wavelength independent; since the 
Fourier transform occurs at the plane where the source is imaged the above 
statement, ri =0. implies that the Fourier transform is also achromatic. This result 
supports the analysis and the analysis itself provides a more general form for 
achromatizing the system for Fraunhofer or Fresnel holography. The next step is to 
describe an optical system which allows the interference of waves of arbitrary 
curvature rather than being restricted to a plane wave and a spherical wave. 

A Three Lens System 

The three lens system of Fig. 4 is capable of forming zone plate holograms 
with beams of arbitrary curvature. It is essentially identical to the four lens system 
without the fourth achromatic lens. L 4 . This is a result of the fact that the final 
achromatic lens serves only to image from one plane to another plane in the optical 
system. Thus, if the light emerging from lens L 4 forms, at some plane, an 
achromatic interference pattern, this plane can be regarded as the image, formed by 
L 4 , of an object oiane elsewhere, in which the two object beams also combine 
achromatica". The '^nn^TcaNon will be different, thus the zone plate recovered a; 
this other ptaoe .v : be o !feront ! y scaled. Also, in general, the two interfering waves 
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The analysis of the three lens system proceeds in the same manner except 
that now both beams have curvature. Expressions for the radius of curvature of the 
interfering wavefronts must be obtained for both beams. Proceeding as before, 
finding the source image with the point source as the object located at infinity, one 
finds the radii of curvature, with reference to Fig. 4, to be 

r 1 = f 3 f 2 2/ ( f 2 2 - f 1 f 3) - d (15) 

r 2 = d (16). 

The plane located by d is the interference plane to be determined by the 
requirements for achrornaticity. Applying these reqirements one first needs the sum 
of the inverses of the radii of curvature which upon substituting wavelength 
dependence becomes 

1 /r^ +-1 / r 2 d 0 1 11 v ' ^30 + (AA/Aq)(1 * 2 d/f 3 g - d A /\ / f 3 0 q ^ (1^) 

We ignore the wavelength squared term, assuming a narrowband approximation, 
yielding the result 

1/ ri + 1/r 2 - (■■ /d - Q >[1 • f 10 d.f 2 2 - d/f 30 f (.V\/A 0 )(1 - 2d/f 30 )] (18) 

The sum is proportional to wavelength only when the plane located by d is at 
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d = f3Q/2. This is the plane where the HOE will be formed. Again the fringe pattern 
will only be approximately achromatic and the recording will be limited to some 
bandwidth. Discussions of these limitations will be found in a later section. 

It is necessary to make a few restrictions on the optical system in order to 
insure that we interfere a converging with a diverging wave. These restrictions are 
taken care of by ensuring that both r-j and r 2 are greater than zero. This is done by 
proper choice of the optical system parameters f-j, f 2 , and fg. With this optical 
system one is able to construct holographic optical elements with a wide range of 
conjugate focal planes and a high signal to noise ratio is obtained by use of a broad 
spectral bandwidth. 

Lim itation s Of Poly chromatic Light 

In designing HOEs in broad spectrum light one must keep in mind the 
consequences of increasing the spectral bandwidth. The space bandwidth product 
(SB) of the HOE is decreased significantly. However, the resulting SB product is 
adequate for most purposes. Using the equation 

( R m 2 + ri 2 ) 1/2 + ( Rm 2 + r2 2 ) 1/2 . fr r2 = (19) , 

where r-| and r 2 are functions of wavelength as in Eq. (15) and (16) and R m is the 
radius of the m th zone, a SB product can be determined. Using numerical 
techniques, we can solve the above equation for R m . For this example we used 
flQ=-15cm., f2=30cm., andf3 0 =15cm. The graphs in Fig. 5 show how R m varies 
with wavelength for each integer m. If we allow for a 1% error in the periodicity of 
R m , that is if 

R m ,max(M * Rm.rrW''-) < [ R m<V * R m-lM(- 01 ) (20) 

we can write the SB product as 

SB = R m (R m " R m-l) (2D- 

With a bandwidth A'. =400 A and \ Q =6000 A we get an SB product of 643. A typical 
polychromatic interference pattern of this bandwidth has a SB product of 







approximately 50 (given by SB = 4\o/A\). The compensated optical system is a 
factor of ten better with an F# of 4.3. There is much to gain in terms of signal to noise 
ratio when using polychromatic illumination. The penalty, in terms of SB product, is 
not too severe. 
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Source-Image Distortion Description of 
Broad Source Fringe Formation 






Inl ro ducticm 

Interferometry with non-coherent light is of great interest to researchers in the 
field of optics 1 ' 29 . The reason for this are twofold; the first is a vastly increased 
signal to noise ratio over use of coherent light and the second is to be able to operate 
interferometers in real world situations (i.e. under ambient illumination). There are 
three classes of interferometers which use non-coherent or partially coherent light. 
First, there are interferometers which operate in polychromatic point source 
illumination 1 '®. Second, are interferometers which operate in illumination of 
reduced spatial coherence 7 ' 17 . And third, there are interferometers which operate 
under illumination that lacks both spatial and temporal coherence 1 ®" 29 . This paper 
discusses the second class of interferometer, that is an interferometer that operates 
under spatially incoherent (or spatially extended source) illumination. 

The interferometers of concern here split the amplitude of the incident light 
into two beams that arrive at the fringe formation plane (output plane) from different 
directions. As seen from the output plane, there are two apparent sources, each an 
image of the actual source. For an extended source, all points on the actual source 
are incoherent with each other, a property that carries over to the two source images. 
However, each point source on the one image has a point on the other image with 
which it is coherent; these two coherent points are each an image of the same actual 
source point. These corresponding point pairs form interference fringes. If extended 
source fringes are to be formed, the coherent pairs should all form identical fringe 
patterns. 

If two such point-pair-coherent sources impinge from different directions, the 
various coherent point pairs produce fringes that generally are of somewhat different 
spatial frequency and may have phase differences that produce fringe patterns that 
are displaced or out of phase with each other. For example, the coherent pairs (a,a 
and b,b') of Fig. 1 may produce fringes of the same spatial frequency, but the fringes 
due to one pair might be shifted relative to those from the other. Also, even if the 
angular separation of the point pairs is the same for each set, the fringes will 
generally not be of equal spacing. For example, suppose the pair aa' produces 
beams impinging symmetrically on the fringe plane P, whereas the point pair bb 
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does not. The fringe spacing will be less in latter the case. This becomes obvious if 
we consider a plane P' perpendicular to the bisector of the beams from b and b'. In 
this plane the fringe spacing will be d, but the fringe pattern in the plane P will be a 
projection, d • cos(AG), where AB is the angle between P and P'. It is evident, then, 
that to have perfect fringe contrast with an extended source, the two apparent-source 
images must be distorted in some appropriate way, so that point pairs from off-axis 
source points have angular separations that increase as a function of the off-axis 
position of the actual source point. 

Our aim is to characterize the required source-image distortion, and to show 
that certain interferometers achieve this distortion exactly, whereas some achieve it 
only approximately, in varying degrees. Conventional interferometers, such as the 
Mach-Zehnder or Michelson, formed from mirrors and beam splitters have no 
distortion at all in their source images, which results in very poor extended source 
source fringe formation capability. We consider symmetrical and asymmetrical 
grating interferometers, as well as a hybrid formed from gratings or mirrors with 
prisms. 

We find that the source distortion factor can be an indicator of the inherent 
extended source fringe forming capability of the interferometer. 

Interferometers that fail to produce the ideal source distortion factor may be 
operated under conditions of reduced spatial incoherence, particularly if a source 
structure factor can be found that optimizes their fringe forming capability. An 
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example of this is the Mach-Zehnder interferometer, whose extended source 
capability was discussed by Bennett 5 and Swanson 26 ; the latter found that the 
Mach-Zehnder has an optimum fringe forming capability when illuminated by a 
continuous ring of point sources that lie at a constant radius from the optic axis. The 
concern here is not with the optimum source structure for a given interferometer, 
rather this paper poses the question: What two source image structures, that have 
point by point coherence between each source, create a perfect parallel fringe 
pattern? It is clear from Bennett's paper and Swanson's thesis that mirrors and beam 
splitters, that produce two exact replicas of the original source at specific locations 
with reference to some interference plane P, do not produce a spatially invariant 
fringe pattern. The ideal optical system that is needed will split the light into two 
beams and then recombine them at some plane to form a set of parallel fringes for a 
spatially extended source of any size. 

I deal Exten ded Sou rc e Struc ture 

The analysis that follows assumes that the illuminating source is located at 
infinity, so that each source point can be considered to be mapped into a plane 
wave of infinite extent. It is a well known tact in coherence thoery that if an infinite 
field is completely incoherent in one plane, it is completely incoherent in every other 
plane. Thus, to the extent that the source can be approximated as infinite, there is 
no loss of generality in assuming a source at infinity, even if the actual source is not. 

Consider two sources in space as shown in Fig. 1. The sources are assumed 
to be located at the back focal plane of a collimating lens with the center of each 
source located symmetrically about the central axis, as seen in Fig. 1. This 
arrangement is to be analyzed with two coordinate systems. Two sets of plane 
waves, each set emanating from a different point on the extended source, have 
been split by some optical system and are brought together at some plane P. These 
plane waves create fringe patterns described by the equations 

IA = 2 + 2cos[4;rxsin(G A )/\ + P A ] (1) 

Ig = 2 * 2cos[47tx'sin(Rg)/A. + Pg] (2) 

G A and Gg are assumed to be different since it is apparent that if they are equal as 
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in the Mach-Zehnder interferometer, the fringe pattern is space variant. The 
consequnce of this is a distortion in the angular location of each source point. This 
distortion is needed in order to obtain a spatially invariant fringe pattern. The source 
point associated with Eq. 1 is assumed to be the zero spatial frequency component of 
the source, i.e. the source point on axis. This simplifies the analysis and is still 
completely general since the source point associated with Eq. 2 could be any source 
point. 

The two fringe patterns must have the same spatial frequency in the x 
coordinate system. The relationship between these two fringe patterns is found by 
rotating the x‘ coordinate system to the x. The result of this gives x' in terms of x as 
x' = xcos(AR) (3), 

where A13 is the angle between the two coordinate systems. From this the 
relationship between the spatial frequencies becomes 

sin(l3 A )cos(Al3) = sin((3g) (4). 

If this relationship is satisfied the spatial frequency of each fringe pattern associated 
with a particular source point is equal when rotated to the x coordinate plane. 

The above condition (Eq. 4) is not enough to ensure that the extended source 
fringe pattern is perfect, i.e., all source elements must form identical fringe patterns. 
Not only must all source elements form fringe patterns of the same spatial frequency, 
but the fringes must also be in registry. Thus, the relative phase between each 
fringe pattern must be zero or some integral number times 2k. Upon examining Eqs. 
1 & 2 this relationship can be expressed as 

P A - = m = 0, ±1, ±2, ±3 .... (5). 

Where P A and Pg are phase constants associated with the pathlength from the 
source to the output plane P. This condition says that the difference in optical 
pathlengths for any source point pair must be zero or an integral number of 
wavelengths. 

If an optical system is able to split and recombine light in accordance with 
Eqs. 4 & 5 it will produce a fringe pattern when illuminated by a spatially incoherent 
source of infinite extent. The distortion needed to produce this fringe pattern is 




evident from Eq. 3. The angular separation between each plane wave pair, 20^ and 
2 0 g, must be different in order to produce coincident fringe patterns in the x 
coordinate system. Typical reflective interferometers, such as the Mach-Zehnder. 
will not produce any distortion in the source and will be limited to smaller sources. 
Refractive optical elements such as a Fresnel biprism can be used to recombine light 
in a manner which distorts the source. Similarly, grating interferometers can do the 
same. Thus, these interferometers have the potential for forming extended source 
fringes, whereas the Mach-Zehnder type, which uses mirrors to achieve their result, 
have considerable limitations in their extended source fringe formation capability 
The Grating Interferometer 

An optical system that will distort the projected source angle can be realized 
with the grating interferometer of Fig. 2. This fringe forming system bends light 
depending on the sine of the angle of incidence and it is this property that causes 
distortion when the sources are projected back as in Fig. 1. The basic grating 
equation is 

sin(Bj) + = M (6). 



Fig 2 Gratings Gi and G2 have spatial frequency fi and ‘>2 respectively. 


The angles in Eq. 6 are seen in Fig. 3, f is the spatial frequency of the grating arid 
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the wavelength of the light being used. 



Fig 3 

The angles that describe the projected source are derived using this basic equation. 
The resultng angular relationships, with respect to Fig. 2, are 


13 1 = sin' 1 [\(f 2 -f 1 ) + sin(l3 0 )] 

(7) 

B 2 = sin- 1 [\(f 2 -f 1 )-sin((3 0 )] 

(8). 

These angles describe plane wave pairs due to a single source point. This source 
point is described by R 0 , its actual angular location, and varies from zero to tninety 

degrees. 


Adapting these last equations two equations to Eq. 4 is an exercise in 

geometry. The results, given below, are easily inserted into Eq. 4 to become Eq. 12. 

f3 A = sin- 1 [^(f 2 -f 1 )] 

(9) 

Rg = (R-l + R 2 )/2 

(10) 

AR = (Ri - R 2 )/2 

(11) 

sin[(Ri + R 2 )/2] • cos[(R 1 - R 2 )/2 = \(f 2 -f 1 ) 

(12) 

Using a trigonometric identity Eq. 12 reduces to 


[sin(Ri) + sin(R 2 )]/2 = \(f 2 -fi) 

(13). 

Substituting R-| and R 2 into Eq. 13 results in an identity. 

Thus, the spatial frequency 
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of the fringe pattern produced by each point pair is exactly equal for the grating 
interferometer. 

As previously mentioned, having the same spatial frequency fringe pattern 
for each point pair is not sufficient for obtaining extended source fringes. The optical 
pathlengths of each arm of the interferometer must have the same pathlength 
measured to the optical axis at the output plane P in order that fringe patterns due to 
each source element are in registry. The distances traveled in each arm of the 
interferometer for any source point is calculated by using the grating equation (Eq. 6; 
The distance traveled (with reference to Fig. 2) by the upper arm is given by 

Di = d r '{1 - \K1, - sin(!3 0 )] 2 } 1/2 + d 2 /{1 - [\(f 2 - U) * s.n(R 0 )] 2 ) 1/2 (13) 

and the lower arm by 

D 2 = di/{1 - [At-, + sin(3 0 )] 2 } 1/2 + d 2 /{1 - [\(f 2 - • s.n(l3 0 )] 2 } 1/2 (14) 

Di and D 2 are exactly equal when d^ is equal to d 2 and f 2 is equal to two f-| In this 
symmetric configuration the grating interferometer is capable of producing perfect 
fringes in extended source illumination. 

The source distortion is achieved by the nature of the gratings and is shown 
graphically in Fig. 4. 



Fig, 4 The graph shows the distortion provided by the grating interferometer to oDta," perfe-- 
extended source fringes. Upper source refers to the upper arm of the interferometer and lower "so-m- 
refers to the ! ower arm 
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Possibilities for source distortion with other optical elements, such as the Fresnel 
biprism, suggest that there are other ways of obtaining perfect fringe patterns in 
spatially incoherent light other than by use of the grating interferometer. 
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Generalization of some incoherent spatial 
filtering techniques 


Emmett N. Leith and Daniel K. Angell 


A two-channel incoherent spatial filtering system is generalized by considering two object transparencies, one 
in each of the two channels. Various special cases result, including two that have previously been described 
and others not previously given. The system can be either linear in irradiance i the basic incoherent easel or 
linear in field ithe basic coherent easel even though the illumination is incoherent in either case In 
particular, we show spatial filtering with the object in one channel and the spatial filter in the other channel. 


I. Introduction 

Spatial filtering with spatially incoherent light has 
attracted considerable interest in recent years. Gor- 
litz and Lanze, Lohmann. Rhodes and Stoner 1 have 
considered two-channel incoherent spatial filtering 
systems, in which the light from the object is split into 
two parts and separately filtered in each before recom¬ 
bining at the output. If H\ and H* are pupil function 
masks placed in each of the two channels, the system 
transfer function is the autocorrelation function of (//, 
+ //■_<>. The terms H j * //, and H 2 * H> are the 
components that result from conventional one-chan¬ 
nel incoherent spatial filtering, while H j « H> and //j * 
H,, the cross-correlation of the two pupil functions, 
result from the presence of two simultaneous channels. 
Since H\ and H< can be any complex function, the 
resulting cross-correlation will generate any arbitrary 
system transfer function It appears that this proce¬ 
dure leads to incoherent spatial filtering systems that 
have the full versatility associated with the coherent 
counterpart. 

Thus, if an object has amplitude transmittance s and 
irradiance transmittance U| J , the cross-product term 
yields an output |.sj- * h, where h is the point spread 
function corresponding to the transfer function H\(i 
H-. Angell" found significant advantages when the 
interferometer is formed from diffraction gratings. 
Collins' among others has noted that by placing the 
object in one of the beams instead of before the inter- 
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ferometer. the optical system becomes linear in ampli¬ 
tude and preserves phase, even though the light is 
spatially incoherent. He demonstrates the construc¬ 
tion of Fourier transform holograms made in spatially 
incoherent light. 

We present here an analysis that generalizes the 
above results. 

II. Basic System 

We show (Fig. 1) an interferometer formed from 
diffraction gratings Gj and G>. The upper beam is 
produced by the +1 order of grating Gj and is redirect¬ 
ed as the -1 order of the grating G>. Similarly, the 
lower beam is produced as the -1 order of Gj and is 
redirected as the +1 order of Gj. The beams are 
brought together again. Although many variants of 
grating interferometers are possible, we chose foi our 
analysis the highly symmetric arrangement where 
gratings G \ and G< have spatial frequencies of /j and 
2/j, respectively. We show optical spatial filtering 
systems in each beam. Without these systems, the 
beams would simply recombine and form, under spa¬ 
tially incoherent illumination (either monochromatic 
or polychromatic), high-contrast fringes of spatial fre¬ 
quency 2/j. When the optical systems are in place, 
fringes become modulated with information related to 
the masks .$j, sj, H j, and H>- The optical system in 
each beam is a spatial filtering system. The upper 
beam contains the object (or signal) transparency spa 
second mask H j placed at the focal plane of the first 
lens is a conventional spatial filter. A similar svstem is 
placed in the lower beam with transparencies ,<_• and 
H j. The two beams are combined at the plane where 
the spatially filtered images of s, and s ■ form Alter¬ 
natively, the two spatial filtering systems could be one 
system with the two beams passing through different 
portions of the same lenses. For the following analv- 
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sis, we chose the former system, although the latter 
yields similar results. 

Ill Theory 

For simplicity, the analysis is carried out in two 
dimensions, x and z. The other lateral dimension \ is 
dropped; however, the analysis we develop is readily 
extended to they dimension also. 

Let the source he extended sufficiently so that it can 
he regarded as completely spatially incoherent, i.e.. an 
infinite source. One element of the source presents a 
plane wave exp(y2T/,,.t i to the grating (?|, where = 
sim'L'A), with the angle of irradiance and A the 
wavelength. The upper beam can be written as 
expt’-"l/.. + f\)x] after exiting from grating G\. Con¬ 
stant coefficients, irrelevant to the analysis, are 
dropped. The light impinging on G _■ has incurred a 
phase delay exp[-;ir\d(/., + f\)-]\ this is the transfer 
function of free space." (We have thus described the 
light propagation in the frequency domain.) A bulk 
phase delay expt j'lxf^d) has been dropped. The field 
exiting G : is then exp;[2;r (/,, - f\)x - 7rAd(/ 0 + /, )-]. 

The delay term ir.Wi/,. + f\ )- is the same for the 
upper and lower branches of the interferometer and 
may. therefore, be ignored. 

Next we consider the field incident on mask si. The 
ma>k and the following optical system are oriented 
perpendicular to the ray from the center of the source; 
thus the term expl -2irjf ] x) is removed. The field 
emerging from the filter transparency is 

//../, i|N.u> ex, »L>2T./ 1) + />|cbr. ill 

where /■ is a spatial frequency variable, f, = xVkF, 
where x' is the physical coordinate of this plane, and F 
is the lens focal length. At the output, the field is 
•; = I Ayr, expL-J../:+ /,•*!*( 

A similar analysis for the lower beam yields 

- .1’//.■/,. f».i ii ex, v'.-T./„ + ; t i|dx 

* expi 131 

Since the two beams combine at nonzero angles we 
have the partial field due to the source element at Do 
being 

u= u, expi-jgT/jxi + u . expi/M,xI '41 

with irradiance l u\ : . The total irradiance is found by 
integrating over the source: 


The first term is just what would result if only the 
upper beam were present, the second that which re¬ 
sults from the lower beam alone. These terms are of 
no interest here. The third and fourth terms are of 
interest. They combine to give a modulated fringe 
pattern of spatial frequency of 2/i. Since they are 
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complex conjugates, we need examine only one of the 
two. The third term can be written 



Fourier transforming the integrals in parentheses 
yields 

* j.i.*xpi . 

which can be rewritten as 


Integrating on / ; , produces a 6-function 6<o - 3) yield¬ 
ing 

This result has several interesting special cases: 

Case 1. s i = S'j = s. This leads to the result 

or. Fourier transforming, we obtain 

.? I = S<H.. H : >. .:;i 

where S is the Fourier transform of si-'. 

This case can be construed as the one already t reated 
by others, where the object '.s| J is placed before the 
interferometer, and thus is presented to both interfer 
ometer branches. The result is an incoherent spatial 
filtering process, where the object irradiance s - is 
modified by a transfer function that is the autocorrela¬ 
tion function of the pupil masks H : and H 
Case 2. s.< = 1, H; = 6-function. We obtain 


In this case, the system is linear in amplitude, just as 
in the case of coherent illumination. However, the 
illumination here is spatially incoherent. In short, the 
incoherent system is behaving like a coherent system 
analogous to the system described by Collins. We 
obtain an arbitrary transfer function without the need 
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for two-pupil masks. However, the system acts on the 
amplitude rather than the irradiance transmittance: 
this difference could either be advantageous or disad¬ 
vantageous depending on circumstances. In particu¬ 
lar, the process would not be effective for diffuse ob¬ 
jects. just as a coherent spatial filter is ineffective for 
diffuse objects. 

Case 3. s- : = 1, Hi = 6-function. We obtain 

/ = >, ff| i Ki> 

This is an interesting and puzzling result. Coherent 
spatial filtering has been carried out with spatially 
incoherent light, as in Case 2. However, we have two 
optical spatial filtering systems in parallel with the 
object in one system and the spatial filter in the other; 
here the light transmitted by the object does not pass 
through the spatial filter, although the effect is as if it 
did. 

If we take H_> = 1, we have a simple imaging process. 
However, the optical system that images the object to 
the output has a vanishingly small aperture, hence, is 
supposedly incapable of image formation with any res¬ 
olution. Yet the resolution can indeed be very good, 
since the effective aperture for the system is not the 
aperture of the object bearing system but is instead 
the aperture of the other branch of the interferometer. 

This seemingly impossible result can be explained in 
a simple way. We think of the source as made up of 
many points. The point element on-axis forms an 
object Fourier transform centered on the pinhole spa¬ 
tial filter, so that only the zero spatial frequency com¬ 
ponent of the object is passed. Another source ele¬ 
ment forms a displaced Fourier transform so that a 
nonzero spatial frequency component /t, with compo¬ 
nents /it,/;., is passed. In this way, each source ele¬ 
ment passes one spatial frequency component of the 
object. A photographic record is then made on the 
output. The zero order of this record, considered as an 
image plane hologram, will form an image with ex¬ 
tremely poor (or no) resolution; however, in the first 
diffracted order the resolution will be to be far in 
excess of what is possible given the pinhole aperture. 
How good the resolution is depends on the size of the 
source. The more extended (i.e., the more spatially 
incoherent) the source, the better will be the resolu¬ 
tion. 

Case 4. s,. s->, H 1t H> arbitrary. This case, which is 
described by Eq. (9). would seem to add little general¬ 
ity to what is offered by Case 2. We can generate any 
possible transfer function by a single transparency Hi. 
So what can he the advantage of the second HC We 
speculate that we obtain better contrast fringes and 
hence a stronger cross-term signal by having two trans 
parencies. For example, for a binary filter H\ =Oorl. 
we might obtain better contrast by having the second 
transparency and setting H : = H 

IV. Experimental Results 

Two experiments were performed to demonstrate 
the theory presented here. The optical system was 




Fig 1. Image ol straightedge formed with jimfrole »|uti.il niter 
conventional one-channel imaging *v-tr-m 


essentially that shown in Fig. 1, except that a single 
lens was used in each beam instead of a pair O’ was of 
100-cycle/mm spatial frequency; 0. was in fact two 
separate gratings, each of 200-cycle mm spatial fre 
quency The source, formed by shining a diverging 
laser beam onto a rotating ground glass was adjustable 
in size and for the most favorable results subtended an 
angle of 4 2° at the collimator. The image was formed 
at the output plane with ;lx magnification 

The first experiment was performed to demonstrate 
Case 2. Edge sharpening was carried out by placing a 
small circular dot in the frequency plane of beam 1 tthe 
object-bearing beam I. and a small pinhole was placed 
in the other. The object was a straightedge. The 
output was recorded as an image-plane hologram and 
read out with a monochromatic broad light source 
sufficiently small that the diffracted orders could be 
separated. The experimental results iFig. 2) show 
prominent edge sharpening, thus verifying the conclu¬ 
sions of Case 2. 

The second experiment was to demonstrate Case 8 
We show improved imaging with the introduction ol a 
second beam, i.e., //_ = ol/,,/.) and s, = 1 The pinhole 
for the experiment was placed in the Fourier transform 
plane of the object beam The second beam had no 
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Fig. 4 Image of straightedge formed with pinhole spatial filter in 
object beam and no restrictive aperture in other beam. 


spatial filter; hence the aperture was effectively the 
image of the source, which was many times larger than 
the pinhole. Figure 3 shows the image of the straight¬ 
edge formed in the conventional way, with only one 
channel. As expected, the image is quite poor because 
of the small pinhole aperture. Figure 4 shows results 
obtained when the other branch of the interferometer 
was in operation, thereby producing a fringe pattern. 
This, on recording, was treated as an image-plane holo¬ 
gram and was read out with a somewhat extended 
nonochromatic light source. The image, forming in 
the first diffracted order, is seen to be considerably 
sharper than in Fig. 3, showing that the effective aper¬ 
ture is far larger than the pinhole, thus verifying the 
conclusions of Case 3. 

This work was sponsored by the Air Force Office of 
Scientific Research grant AFOSR-85-0012 and by the 
National Science Foundation grant NSF-G-ECS- 
8212472. 
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6. Superresolution bv Incoherent to Coherent Conversion 

The previous section discussed a variety of situations in spatial filtering with 
incoherent light. We showed that if the object is placed in only one branch of the 
interferometer (as in Fig. 1), the system acquired a number of interesting properties. 
We are always interested in the cross-product term, recorded on a spatial carrier, as 
in holography. With a signal s-| and a spatial filter H-| in one branch, and a signal S 2 
and a spatial filter H 2 in the other branch (as in Fig. 1), the cross-product term, which 
is read out as the 1st diffracted order of the recorded output, has the form 
S 1 S 2* * h 1 h 2* 0) 
where the * denotes convolution. 



If we take Si = S 2 = s. we obtain immediately the result 
|s| 2 * (h 1 h 2 *) (2) 

which is equivalent to the well-known case of having the object upstream from the 
interferometer. The system is linear in intensity, and operates with object irradiance 
jsl 2 . The transfer function H is the Fourier transform of the point spread function 
h 1 h 2 *. or 

H = H 1 ©H 2 (3) 


where ©means cross-correlation. 







Other cases not as yet explored, are also possible. An interesting possibility is to 
have S 2 = 1 and h 2 = 1 whereupon we have an output 
St * hi (4) 

This is mathematically the case of coherent spatial filtering, even though it has been 
carried out with spatially incoherent light. This interesting possibility is one we wish 
to explore. We can let Hi be any filter that has been made for coherent illumination. 
It can, for example, be a holographic spatial matched filter, or a complex filter image 
restoration, or deblurring. We expect to obtain all of the advantages of coherent 
spatial filtering, along with the low noise excellent image quality associated with 
incoherent imaging. 

Another variation of this, which is mathematically equivalent, but physically quite 
different, is to have hi = 1. The operation is then 

si * h 2 * (5) 

As before, this is a coherent operation carried out with incoherent light. The object 
Si has been spatially filtered by a filter H 2 that is physically isolated; the light that 
emanates from the object Si does not pass through the mask at all, yet the mask 
exerts the same effect as if it were in the object light path. Here is a new concept in 
spatial filtering, where the signal and filter are essentially isolated. 

Again, suppose s 2 = 1. Let the filter impulse response be binary (0 or 1), with 
hih 2 * = hi = h 2 = h. Is there any detectable difference between having a spatial 
filter H in just one path or both? The mathematics suggests identical situations, but 
of course physically they are quite different. What advantages can arise from having 
H in the other channel or in both channels? Would the SNR be improved? 

The most parodoxical case of all comes by letting s 2 = 1, hi =1 and h 2 the 

desired spread function. Note that hi = 1 means Hi = 6 (a dirac delta function); this 
situation is depicted in figure 2. 








Fig. 2 

If we consider the upper path as an imaging system by itself, it is a system 
incapable of much resolution. If the pinhole is indeed a delta function, i.e., truly a 
point, there is no resolution at all. If the pinhole is of a more practical size, say 10-50 
microns, which is a common pinhole size, then the resolution will be extremely poor. 
This situation is equally true for either coherent or incoherent illumination - the 
resolution is not much different for these two cases. Indeed, within the framework of 
conventional imaging theory, the resolution simply cannot be better than the 
Rayleigh limit imposed by the aperture size. 

Yet in this case theoreticoal considerations predict resolution equal to what one 
would get from an imaging system without the pinhole constriction. We thus obtain 
resolution far in excess of the classical limit - in short, superresolution. 

One might ask, how can the resolution be that of a normal imaging system? The 
information needed to produce the resolution is simply not there - such information is 
not transmitted by the aperture. The aperture of the parallel system is adequate, but 
then, no object information passes through this aperture. 

Techniques exist for superresolution. One common method is analytic 
continuation. Such techniques require nearly noise-free systems, extremely 
accurate measurements, and considerable computation. They can often be used to 
go a factor two beyond the classical limit, but seldom any further. 

Are we dealing here with a similar superresolution system that would have the 
same limitations, or is it something quite different? We suggest the latter, and our 
initial experimental results support this view. Certainly, the exacting measurement 
and computational aspects of superresolution are absent here. 

What would one expect the resolution to be in a practical system? In practice, the 
source will not be infinite, but of a size of the order of a centimeter. This source is 
imaged in beam two (Fig. 2). Now, suppose we place an iris in beam two at this 
source image plane. Let the iris be initially small, and we proceed to open it. As we 
do so the resolution is expected to improve. When the iris is the size of the source 







image, further enlargement lets through no additional light. Thus the effective 
aperture of the system is the size of the source image. Or put equivalently, the 
effective f-number of the imaging system depends upon the angle the source 
substends at the collimator. The larger the source, the larger the effective aperture, 
and the better the resolution. 

Carrying this viewpoint further leads to the conclusion that this system is 
essentially different from other superresolution techniques. We show that the gain in 
resolution does come at a price, but at a fairly modest one. 

An incoherent (or partially coherent) system can be considered as the 
superposition of many coherent ones; each source element gives rise to a separate 
channel with its own independent noise. If the optical system remains the same, the 
number of resolvable elements in the image is the number of effective coherent 
sources, and hence the number of independent channels. One would expec* the 
incoherent system to have a SNR improvement of the N 1 '^ over the coherent case, 
where N is the number of channels. In practice, the improvement is enormously less 
than this, since the noise produced by adjacent or nearby channels would not be 
independent, particularly for noise of components of low spatial frequency. 

The above optical system in effect converts the incoherent system into a coherent 
one in the following way. The center element of the source forms, at the pinhole 
aperture of channel 1 , a Fourier transform of the object s. The pinhole allows only 
the 0 spatial frequency to pass. Similarly an off-axis element projects a displs or 
off-axis Fourier transform to the pinhole, and a nonzero spatial frequency .asses. 
Thus, the number of spatial frequency components that pass is the number of 
resolvable source elements as seen from the Fourier transform plane. For an 
incoherent system, these contributions would combine incoherently and would loose 
their identity as different spatial frequency components of the object. But in the 
above case, they combine coherently, even though the source elements are 
completely incoherent. The effect is that of a source the size of pinhole, and an 
imaging aperture the size of the source. 

This tradeoff can often be highly favorable, since experience shows that 
moderate reduction of redundancy due to incoherency generally leads to no 
observed degradation of the image. The number of redundancy channels in the 
output system can be of the order of 10® and reducing this to an order 10^ typically 
does not make much difference, either in practice or in theory In practice the 
difference is far less than in theory, since closely adjacent channels carry only 
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partially uncorrelated noise, ie.e, low spatial-frequency noise components would be 
almost identical for adjacent channels. 

This appears to be a basically new concept in superresolution: superresolution 
by incoherent to coherent conversion. The discovery of this concept was rather 
serendipitous, and we have since pondered on what its implications may be, and on 
what usefulness it can have. We raise several pertinent questions: 

1. Can it be used to overcome aberrations? For example, the imaging system 
carrying the information ,s operated at high f-number (a consequence of the 
pinhole) and therefore will have low aberration. The reference imaging system 
operates at low f-number, but with a very narrow field. Hence, the overall 
imaging process may be better than the comoponent lenses would ordinarily 
allow. 

2. Can the process be achromatized, so that it can be carried out in polychromatic 
light with the possibility of tradeoff of temporal incoherency for resolution? 

3. How can the basic system be modified and extended so as to lead to practical 
applications? 

4. Can we use the system in a situation where high light levels, used in one branch, 
can be spatially filtered with low light levels in the other branch? 

7. Journal Articles Prepared or Published This Year 

1. S. Leon and E. Leith, "Optical Processing and Holography with Reduced 
Coherence Polychromatic Point Source Illumination," Appl. Opt. 24, 3638 
(1985). 

2. E. Leith and D. Angell, "Generalization of Some Incoherent Spatial Filtering 
Techniques," Appl. Opt. 25, 499 (1985). 

Some papers almost ready for submission 

3. S. Leon and E. Leith, "Generalization of the Two-Zoneplate Method of 
Achromatized HOE Formation," (Cpt. 3 of this report). Will shortly be 
submitted to Applied Optics. 

4. S. Leon and E. Leith, "Source-Image Distortion Description of Broad Source 
Fringe Formation." (Cpt. 4 of this report). Will shortly be submitted to Applied 
Optics. 
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5. E. Leith and D. Angell, "Superresolution by Incoherent to Coherent 
Conversion," (Cpt. 6 of this report). Will shortly be submitted probably to j. 
Opt. Soc. Am. 
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E Leith (PI) 

S Leon (Graduate Student) 

D. Angell (Graduate Student) 


Of our various accomplishments, one stands out as being of special interest. This 
is our discovery of a new type of imaging system, in which the optical transfer 
function of the imaging process is controlled by the optical transmission 
characteristics of a parallel channel, through which light from the obiect does not go 
A dramatic example of this is obtained by using a very small aperture in the imaging 
channel, and a larger aperture in the other channel. Resolution far beyond the 
diffraction limit is obtained. We believe this to be an astonishing result, wholly new. 
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